AO-A042  644 


UNCLASS IF ICO 


ROCHESTER  UNIV  N T DEPT  OF  STATISTICS  F/t  l*/l 

testing  significance  of  a mean  vector  - A POSSIBLE  ALTERNATIVE  — ETC(U) 
1977  G S MUOHOLKAR*  P SUBBAIAH  AFOSR-3360-77 

AFOSR-TR-77-0856  NL 


AOA042644 


AFOSR-TR-  7 7 - 0 8 5 6 

TESTING  SIGNIFICANCE  OF  A MEAN  VECTOR - 
A POSSIBLE  ALTERNATIVE  TO  HOTELLING'S  T“  * 


by 

Govind  S.  Mudholkar  and  Perla  Subbaiah 
University  of  Rochester  and  Oakland  Univers 


D D C 


Abstract 

In  problems  involving  multivariate  measurements  experimental 
considerations  often  indicate  grouping  of  variables  into  subsets  ordered 
according  to  their  importance.  In  such  situations,  the  problems  such  as 
comparison  of  two  mean  vectors  and  profile  analysis  may  be  treated  by 
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Hotelling's  T“-test  adapted  along  the  lines  of  the  step-wise  procedure 
of  J.  Roy  (1958),  or  the  well  known  test  for  additional  information  due  to 
Rao  (1948).  In  this  paper  we  study  a modification  of  the  step-wise  pro- 
cedure obtained  by  combining  the  component  tests.  The  exact  Bahadur 
slopes  of  resulting  procedures  are  computed  and  it  is  shown  that  the  pro- 
cedure based  upon  Fisher's  combination  method  is  asymptotically  equivalent 


to  Hotelling's  T".  A Monte  Carlo  study  suggests  that  even  in  small  sam- 
ples the  power  functions  of  the  new  method  and  Hotelling's  T"-test  are 
practically  equivalent. 
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1-  INTRODUCTION.  Hotelling's  T -test  which  involves  all  the  variables 
symmetrically  is  the  most  common  method  of  testing  multivariate  hypotheses 
such  as  equality  of  two  mean  vectors  or  similarity  of  the  profiles  of  two 
groups.  This  may  not  be  appropriate  if  the  variables  are  of  unequal  impor- 
tance as  in  manv  biological  experiments  where  the  measurements  are  often 
associated  with  biological  processes  and  are  gouped  into  subsets  which  can 
be  ordered  according  to  their  biological  significance.  A grouping  of 
variables  in  two  subsets  occurs  naturally  in  most  investigations,  the  first 
group  comprising  the  variables  of  primary  interest  and  the  second  being  the 
group  of  less  relevant  variables  obtainable  at  little  additional  cost.  A 
common  procedure  for  testing  the  hypothesis  in  such  situation  is  the  step- 
wise procedure  related  to  the  well  known  test  for  additional  information 
(J.  Roy  (1958),  Rao  (1948)).  In  this  procedure  the  hypothesis  on  the  means 
of  the  variables  in  the  first  group  is  tested  by  the  usual  T“-test,  but  the 
hypotheses  on  the  means  of  the  subsequent  groups  are  tested  by  T“-tests  in 
which  the  previous  groups  are  regarded  as  concomitants.  An  advantage  of 
this  procedure  consists  in  the  independence  of  the  T^-statistics  under  the 
overall  hypothesis  which  makes  the  control  of  type  I error  manageable. 
Customarily  the  individual  tests  are  conducted  at  suitable  levels  and  the 
overall  hypothesis  is  rejected  if  at  least  one  of  the  component  tests  is 
significant.  Alternatively,  one  may  summarize  the  step-wise  procedure  by 
reporting  the  P-values  of  the  component  tests  and  combining  (e.g.  Oosterhoff 
(1969))  them  to  obtain  the  overall  significance  probability. 

In  Section  2 we  present  variations  of  the  step-wise  procedure  related 
to  several  combination  methods  and  discuss  an  underlying  invariance  structure 
leading  to  a canonical  form.  In  Section  3 the  exact  slopes  of  these  varia- 
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tions  are  obtained  for  computing  Bahadur's  ARE's.  In  Section  4,  a Monte 

Carlo  study  and  its  conclusions  are  presented.  The  modified  step-wise 

procedure  based  upon  Fisher's  combination  method  is  seen  to  be  an 

2 

asymptotic  equivalent  of  Hotelling's  T . This  observation  is  supported 
by  the  simulation  study. 


2.  SOME  MODIFICATIONS  OF  THE  STEP-WISE  PROCEDURE.  Let  X,  .X., , . . . ,X_ 
— ^n 

be  n independent  observations  on  a random  vector  X having  a p variate 
normal  distribution  with  mean  y and  covariance  matri.x  I,  and  consider 
the  problem  of  testing  H^:  u = 0.  If  ^ = i^X^,  S = J^(X^  - X)(X^  - X)'  , 
then  Hotelling's  test  rejects  H^  for  large  values  of  T"  = n(n  - 1)X'S  X, 
where  under  H^,  (n  - p)T^/[p(n  - 1)]  = (n  - p)nX'S  ^X/p  is  distributed 
an  F-variable  with  p and  (n  - p)  degrees  of  freedom.  In  case  the  varia- 

T 2 . . 

bles  have  an  a priori  order,  and  TT  denote  the  Hotelling's  T“-statistics 
for  the  first  i variates,  i = l,2,...,p,  then  the  step-down  procedure  for 
MANOVA  specialized  to  this  problem  consists  of  p tests  based  upon  statis- 
tics 


F.  = (n  - DtT^  - T^_^]/[(n  - 1)  ^ T^_^] 


(2.1) 


i = 1,2 p,  and  rejects  H^  if  any  of  the  component  tests  is  signifi- 

cant. The  type  I error  control  of  the  procedure  uses  the  fact  that  under 
Hq  , F^  are  independently  distributed  as  F(l,n-i)  variates,  i = l,...,p. 

The  logic  of  the  step-down  procedure  extends  as  well  to  the  more 
general  and  common  situation  where  the  variables  are  grouped  into  K 
subsets  and  the  subsets  are  ordered  according  to  their  importance.  Let 


the  number  of  variates  in  the 


. th 
1 


subset  be  p. 


and 


j-1 
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q=  Zp.  =p,  i=l,...,k.  The  random  vector  X and  the  parameters 

^ i=l  ^ 

of  its  distribution  may  then  be  partitioned  as 
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(2.2) 


^i  - ^il'^l 


6i2U- 


- ®i  1^-r 


It  may  be  noted 


that  u = 0 is  equivalent  to  the  conjunction  of  6^  = 0,  i.e., 

Hq  = n ^oi’  ^oi»  which  may  be  referred  to  as  the  hypothesis  con- 

ceming  the  "additional  information"  provided  by  the  i’^  subset,  can  be 
tested  using 

Fi  = (n  - qj)|T^  . T )/([(„.  1)  . ]p}  . 

^1  ^1-1  ^i-1  ^ 


(2.3) 


If  H = n H is  true  then  the  F.'s  are  independently  distributed 
i=l  ^ 

as  F - variables  with  (p^,  n - q^)  d.f..  The  step-wise  procedure,  in 

this  case  consists  of  K tests  with  critical  regions  F^  ^ F(Pj^*  n - q^;  a^). 


where  F(p/,  n - q^;  a^)  denotes  the  (1  - Oj^)  100^*^  percentile  of 

F(p.,  n - q.)  variate,  i = l,...,k,  and  rejects  H at  level  a = 

i k 1 0 

1 - 11(1  - o. ) if  at  least  one  of  the  tests  is  significant. 

i*l  ^ 


the 


An  alternative  to  comparing  F^  with  F(p^,  n - q^;  o^) 


and 
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labelling  it  as  significant  or  insignificant  is  to  report  the  P-value, 

= Pr(  F(P^  , n - q^)  | associated  with  it  as  the  summary  of  the 

test  of  In  addition  to  avoiding  the  problem  of  having  to  select  the 

levels  of  the  component  tests,  this  approach  permits  an  assessment  of  the 
overall  significance  of  the  data  by  combining  the  P-values,  because  under 
Hq  the  P-values  are  independently  uniformly  distributed.  A combination 
statistic  is  usually  a simple  function  iJ/(Pj^ , . . . ,Pj^)  of  the  P-values  with  a 
.simple  null  distribution  (Oosterhoff  (1969)).  In  the  sequel  we  investigate 
two  such  statistics  );p  = --  A log  P^  and  = min  P^  , where  large 
values  of  *'p  and  small  values  of  indicate  significance.  Under 

'J)_  is  distributed  as  x*"  with  2k  degrees  of  freedom  and  the  distribution 

r 

of  '1/.J.  is  given  by  Pr(iii.j.  ^ C j Hq)  = 1 - (1  - c)  . A summary  of  the  modified 
step-down  procedure  consists  of  the  k P-values  together  with  the  P-value 
of  the  combination  statistic. 

Now  we  present  an  invariance  reduction  of  the  problem  leading  to  a 
canonical  form  which  permits  investigation  of  the  properties  of  the  modified 
step-down  procedures.  It  is  well  known  that  Hotelling's  T“  is  a maximal 
invariant  under  the  group  G of  nonsingular  transformation  of  the  p 

variables,  and  is  UMP  invariant  for  testing  the  hypothesis  u=  0 vs. 

2 

u 0.  The  power  of  the  T -test  involves  only  the  noncentrality 
parameter  y'l  ^u.  If  all  the  variables  can  be  arranged  in  a strictly  de- 
creasing  order  of  importance,  it  is  known  that  the  step-down  statistics  F^, 
i » l,...,p  given  in  (2.P  are  maximal  invariants  under  the  group  of  lower 
triangular  transformation  of  the  p variables  (Subbaiah  and  Mudholkar  (1969)). 
In  the  following  theorem  this  invariance  reduction  is  extended  to  the  case 
of  block-structure. 
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Theorem  2.1.  Let 


be  the  group  of  nonsingular  lower  block  triangular 


matrices  L = (L..),  where  L.  . is  of  order  p.  xp.  and  L.  . = 0 for 

i<  j,  j = 1,2 k.  Then  the  problem  of  testing  Hq.-  u = 0 vs. 

u 0 is  invariant  under  transformation  X -►  LX,  S -►  LSL' , and  the 
step-down  statistics  defined  in  (2.5)  are  ma.ximal  invariants. 

Proof:  The  invariance  of  F,,...,F,  follows  trivially  from  the 

A K 


invariance  of  T“  , i = 1 , . . . ,k.  In  order  to  see  the  maximal  invariance, 

_ 

suppose  that  X , S and  Y , S , the  sufficient  statistics  from  two  data 

« ^ X ^ y 

sets,  give  rise  to  the  same  step-down  statistics.  It  can  be  shown  easily 

that  there  exists  a lower  block  triangular  matrix  L such  that  Y = LX 

and  S = LSL'.  Suppose  L and  L are  lower  block  triangular  matrices 

such  that  S = L L'  , S = L L'  , and  U = L X = (u,  • . • u,  ) ' , V = L'^Y  = 
-X  ,x.  .y  ,y.y  , ,x,  1 k , ,y  , 

(Vj ^k^ ' ■ equality  of  the  step-down  statistics  from  both  data 

sets  indicate  that  u!u.  = v.'v.,  which  implies  that  there  exists  an  ortho- 

gonal  matrix  such  that  M^Cp^xp^)  such  that  v^  = M^u^,  i = l,...,k.  By 

taking  L = L^ML^^,  where  M is  a block  diagonal  matrix,  with  M^,...,N1^ 

as  diagonal  blocks,  it  can  be  shown  that  Y = LX  and  S = LS  L'.  Hence 

,y  - -X  - 

the  theorem. 


Theorem  2.2.  The  lower  function  of  any  invariant  test  of  vs. 

depends  upon  parameters  5^  * ^i^i’  ^ * l,...,k,  where  n = = 

(n^ I • • • ,nj^) ' and  B is  a lower  block  triangular  matrix  such  that  I » BB'. 

Proof:  The  theorem  follows  from  Theorem  2.1,  replacing  X by  u and 

S by  2,  and  noting  that  i5^'s  are  maximal  invariants  in  the  parametric 
space  under  the  induced  group  of  transformation  (Lehmann  (1959)). 
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3.  BAHADUR  ARE'S  OF  THE  MODIFIED  STEP-DOWN  PROCEDURES.  Let  T be 
n 

a statistic  used  for  testing  a null  hypothesis  H^:  9e  ®0  vs.  an  alterna- 
tive Hj^:  9 e ^ , where  large  values  of  T^  indicate  significance.  Then 

the  rate  of  decrease  to  :ero  of  the  P-value  L (t  ) = Pr(T  > t 1 H.) , 

n n n — n ' 0 

evaluated  at  t = T , as  n increases  is  taken  as  a measure  of  efficiency 
n n 

of  the  test.  The  following  Definition  3.1  and  Theorem  3.2  summarize  the 
concept  of  e.xact  slope  and  a useful  method  for  its  computation. 

Definition  3.1.  The  exact  slope  cC9)  of  {T^}  is  given  by 


- ic  (9)  = lim  n'^log  , 
n-*« 


(3.1) 


providing  that  the  (a.s.)  limit  exists. 

Theorem  3.2.  (Bahadur  (1971,  p.  27)).  Suppose  that 

limi,  nS"  =b(9),  a.s. 

for  each  9 c where  -®<  b(9)  < ”,  and  that 


(5.2) 


lim  n*^  log[l  - F^(i4i’t)]  = - f(t)  , 


(3.3) 


for  each  t in  an  open  interval  I,  where  f is  a continuous  function 
on  I.  and  {b(9);  9 e H } <C  I . Then  (3.1)  holds  with  c(9)  = 2 f(b(9)) 


for  each 


9 e 


Now  consider  the  problem  of  testing  the  multivariate  hypothesis 

H^:  u = 0 vs.  H^:  u 0 described  in  Section  2.  The  following  theorem 

2 

gives  the  exact  slope  of  the  T -test  for  this  problem. 

2 

Theorem  3.3.  The  exact  slope  of  Hotelling's  T test  is  given  by 


* log(l  2 n!  . 

i«l'  ' 


(3.4) 
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Proof:  If  we  denote  = [(n-p)nX'S  ^X/p]^^,  then 

_S-S.  1 

T /»n  ? (u'  r u/p)'^  = ( I n.n./p)  ^ and  lim  n log  [1  - F (/n  t) 

i=l'  ' n-HJo 


i=l'  ' n-HJo 

lim  n ^ log  [1  - Pr(n  X S ^X  ^ npt^/(n  - p) ] = - jlog(l  + p t^)  (follows  from 
n-*“ 

Bahadur  (1971,  p.  13).  Hence  the  theorem. 

Now  we  obtain  the  exact  slopes  of  the  modified  step-down  procedure 


based  upon  Fisher's  and  Tippett's  methods. 

Lemma  5.4.  The  exact  slope  of  the  i^^  component  test  is  given  by 


c.(9)  = logfl  + 


1 E n '.  n • 


(3.5) 


Proof:  The  theorem  follows  easily  by  noting  that  for  T = /F.(n), 

' n 1 


lim  n'^'^T^  = {n!n./[(l  > E njnJPi]  )''^  . 

n-«  ' ' j = 


lim  n ^ log  [1  - F (.^t)]  = lim  n ^ log  Pr(F.  ^ nt“) 

n-Hx>  n-«o 


.1  X“(p.)  nt“p. 

= lim  n log  ^ 


X“(n-q.)  (n-q.) 


-|-log(l  ^ t^p.) 


(follows  from  Bahadur  (1971,  p.  13,  equation  5.7)) 


Theorem  3.4.  The  exact  slopes  Cp  , c^  of  Fisher's  method  and  Tippett's 


method  of  conbining  the  step-down  P^'s  are  given  by 


-8- 


k 

'f  ° 1^1  ■ “>s»  • 


= maxc^  = max  log  1 + 


i-1 

1 + E n '•  r| , 


Proof : This  theorem  follows  from  the  results  due  to  Littelland  Folks  (1971) 

A comparison  of  Theorem  3.3  and  Theorem  3.5  shows  that  the  modified 
step-down  procedure  based  upon  Fisher's  method  of  combination  of  tests  is  an 
asymptotic  equivalent  of  Hotelling's  T“  in  the  sense  of  Bahadur  ARE.  On 
the  other  hand,  in  this  sense  the  modified  step-down  procedure  based  upon 
Tippett's  combination  method  is  in  general  less  effective  and  never  more 
effective  than  the  T“-test. 


4.  POWER  FUNCTIONS  OF  THE  MODIFIED  STEP-DOWN  PROCEDURE  BASED  UPON 


A SIMULATION  STUDY.  In  this  section  we  summarize  a simulation  experiment 

conducted  in  order  to  understand  the  moderate-size  sample  behavior 

of  the  modified  step-wise  methods  in  relation  to  Hotelling's  T“.  In  this 

e.xperiment  we  study  the  special  case  pj^  = p^  = . . . = Pj^  = 1 . Our  objective 

is  to  obtain  a relatively  detailed  profile  of  the  power  function  when  p = 2, 

and  an  indication  of  its  general  behavior  in  certain  directions  when  p = 3,4. 

In  view  of  the  invariance  structure  in  the  problem  we  take  without  any 

lose  of  generality  E = I,  in  which  case  the  power  functions  of  the  modified 

step-down  procedures  as  well  as  that  of  Hotelling's  T“-test  depend  only  upon 

2 2 

the  noncentrality  parameters  i =!.-.•• -,p-  For  the  case  p = 2, 

the  power  functions  of  the  methods  are  estimated  over  the  entire  plane 
0'^  , ^2^'  whereas  for  p » 3,4,  the  estimates  are  obtained  for  certain 
directions  only,  namely,  the  equiangular  line  “ ^2  * * '^p  along 
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a coordinate  axis  i.e.,  for  alternatives  (u,0,...,0). 

Monte  Carlo  Experiment:  The  standard  normal  deviates  are  generated 

on  the  IBM  360/365  computer  at  the  University  of  Rochester  using  "McGill 
University  random  number  package"  based  upon  the  technique  of  Marsaglia 
(1961)  for  generating  standard  normal  deviates.  A random  observation  from 
a p-  variate  normal  population  , 1^)  is  obtained  by  drawing  p 

random  observations  from  a standard  univariate  normal  population  and  adding 
.u,  I • • • ,Up  to  them  respectively.  When  p = 2,  the  deviates  are 
generated  for  the  values  of  = 0.0(0. 1)1.6  and  u-,  = 0.0(0. 1)1.9.  IVhen 
p = 3,4,  they  are  obtained  for  (u , ^ , • • • , u)  and  (u ,0, . . . , 0) , for 
u = 0.0(0. 1)1.0. 

For  each  of  the  parameter  values  3000  samples  of  size  n = 20  are  ob- 
tained, and  are  then  used  to  estimate  the  power  functions  of  various  tests 
of  Hq.'  u = 0.  Specifically,  for  each  of  the  samples  we  compute  Hotelling's 
T“  statistic  and  -Jj  and  the  statistics  for  the  modified  step-down 

procedures  related  to  Fisher's  and  Tippett's  combination  methods.  The  IMSL 
routine  MDFD  is  used  for  obtaining  the  P-values  needed  in  the  computation 
of  , itp.  The  values  of  the  statistics  for  each  sample  are  compared  with 
the  corresponding  critical  constants  for  a = .01,  .05  and  .10.  The  power  of 
a procedure  with  a given  y and  a is  estimated  by  the  proportion  p of 
times  Hq  is  rejected  in  the  3000  trials;  the  s.e.  of  the  estimate  being 
(p(l  -p)/3000)^  ^ .009.  The  exact  power  of  Hotelling's  T is  also  computed 
using  the  FORTRA.N’  routine  by  Bargmann  and  Ghosh  (1964)  for  computing  the 
c.d.f.  of  noncentral  F distribution. 

Results:  A selection  of  the  results  of  the  simulation  study  is  given 

in  Tables  1 through  5.  Tables,  1,  2,  3 contain  a relatively  detailed 
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profile  of  the  power  functions  for  p = 2.  In  Tables  4,  5 we  give  the 

empirical  power  functions  of  the  three  tests  and  the  exact  power  of  the  T” 

test  for  p = 3,4,  and  a = .05  corresponding  to  the  two  configurations 

in  the  parametric  space,  vie.,  the  equiangular  configuration  (p,a,...,u) 

and  the  extreme  configuration  (u,0,...0).  The  results  of  the  study  seem 

generally  supportive  of  the  conclusions  drawn  from  Theorems  5.5  and  3.5. 

2 

Specifically,  (1)  the  power  functions  of  the  T -test  and  the  modified 
step-wise  test  based  on  Fisher's  method  appear  indistinguishable  for  the 
cases  simulated,  (2)  the  modified  test  based  upon  Tippett's  method  seems 
to  have  an  advantage  over  the  T^-test  along  the  coordinate  axis.  .-Mong  the 
equiangular  line  the  T“-test  dominates  it. 
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TABLE  1 


EX.ACT  POWER  FUNCTION  OF  HOTELLING'S  T”  TEST  PROCEDURE 


1 = .05,  p = 2 


1.5 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.2 

.995 

.996 

.996 

.999 

.999 

1.00 

1.00 

1.00 

1.0 

. 965 

.966 

.970 

.989 

.998 

1.00 

1.00 

1.00 

0.3 

.845 

.351 

.367 

.943 

.990 

.998 

1.00 

1.00 

0.5 

. 436 

.451 

.495 

.742 

.943 

.989 

.999 

1.00 

0.2 

. 105 

.113 

. 165 

.495 

.867 

.970 

.996 

1.00 

0.1 

. 063 

.076 

.119 

.451 

.851 

.966 

.996 

1.00 

0.0 

.050 

. 063 

.105 

.436 

.345 

.965 

.995 

1.00 

o 

o 

0.  1 

0.2 

0.5 

0.8 

1.0 

1.2 

1.5 

TABLE  2 

POWER  FUNCTION  OF  FISHER'S  COMBINATION  OF  THE  STEP-DOWN  TESTS 
ESTI.MATED  FROM  THE  MONTE  CARLO  EXPERIMENT* 


a = . 05 , p = 2 


1 . 5 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.2 

.994 

.993 

.995 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

.954 

.960 

. 966 

.988 

.999 

1.00 

1.00 

1.00 

0.8 

.324 

.338 

.848 

.947 

.991 

.999 

1.00 

1.00 

0.5 

.416 

.441 

.46" 

.764 

.945 

.989 

.998 

1.00 

0.2 

. 100 

.112 

. 163 

.510 

.864 

.965 

.996 

1.00 

0.1 

.064 

.087 

.115 

.422 

.838 

.964 

.995 

1.00 

0.0 

.050 

. 063 

.116 

.434 

.842 

. 963 

.963 

1.00 

0.0 

0.1 

0.2 

0.5 

0.8 

1.0 

1.2 

1 . 5 

TABLE  3 


POWER  FUNCTION  OF  TIPPETT'S  COMBINATION  OF  THE  STEP-DOWN  TESTS 
ESTIMATED  FROM  THE  MONTE  CARLO  EXPERIMENT* 


a = .05,  p = 2 


1 . 3 

1 .00 

1,00 

1 . 00 

1 . 00 

1 . 00 

1.00 

1.00 

1.00 

1.2 

.99* 

.996 

.994 

.997 

.998 

1.00 

1.00 

1.00 

1.0 

.966 

.9h9 

.964 

.969 

.990 

.998 

1.00 

1.00 

0.3 

. S48 

. 849 

.850 

.890 

.970 

.993 

.999 

1.00 

0.5 

.441 

.446 

.453 

.n~4 

. 9ob 

.981 

.998 

1.00 

0.2 

.101 

. 1 14 

. 159 

. 4*4 

.884 

.969 

.997 

1.00 

0.1 

.06  3 

.073 

. 109 

. 435 

.856 

.974 

.99* 

1.00 

0.0 

.052 

.060 

. 106 

. 457 

.369 

.975 

.998 

1.00 

0.0 

0.  1 

0.2 

0.5 

0.8 

1.0 

1.2 

1.5 

^1 

Each  estimate  is  based  upon  3000  trials. 
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TABLE  4 

POWER  FUNCTIONS  OF  T‘-TEST  AND 
MODIFIED  STEP-DOWN  TESTS* 
a = . 05  , p = 3 


Configuration 

Tests 

u 

0.0 

0.1 

0.2 

0.4 

0.6 

0.8 

1.0 

Fisher 

.050 

.079 

. 195 

. 645 

.953 

.998 

1.00 

Tippett 

.046 

.082 

. 174 

.496 

.844 

.982 

.999 

u 

T“ 

.049 

.078 

. 192 

.629 

.948 

.998 

1.00 

T”-exact 

.049 

.078 

. 192 

.629 

.948 

.998 

1.00 

/ u \ 

Fisher 

.052 

.059 

.095 

.239 

.510 

. 765 

.917 

o] 

Tippett 

.048 

.056 

.090 

.254 

. 559 

.828 

.958 

\o) 

T" 

.050 

.058 

.091 

.234 

.510 

. 772 

.926 

T^-e-xact 

.050 

.060 

.091 

.238 

.498 

.768 

.931 

TABLE  5 

POWER  FUNCTIONS  OF  T“-TEST  AND 
MODIFIED  STEP-DOWN  TESTS* 
a = . 05  , p = 4 


Configuration 

u 

0.0 

0.1 

0.2 

0.4 

0.6 

0.8 

1.0 

Fisher 

.045 

.086 

.01 

.735 

.979 

1.00 

1.00 

Tippett 

.047 

.084 

. 164 

.508 

.837 

.982 

.999 

\“/ 

T^ 

.046 

.035 

.197 

.712 

.973 

.999 

1.00 

u 

T”-exact 

.050 

.082 

.202 

.695 

.973 

1.00 

1.00 

y 

Fisher 

.053 

.064 

.096 

.209 

.436 

.694 

.880 

/ 0 \ 

Tippett 

.048 

.060 

.084 

. 229 

.513 

.779 

.945 

»/ 

T" 

.055 

.061 

.090 

.203 

.427 

.696 

.885 

T'"-exact 

.050 

.058 

.082 

.202 

.428 

.695 

.388 

* Each  estimate  is  based  upon  3000  trials. 
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Additional  Tables 
for 


TESTING  SIGNIFICAiNCE  OF  A ME.AN  VECTOR  - 
A POSSIBLE  ALTERNATIVE  TO  HOTELLING’S  T^  * 

by 

Govind  S.  Mudholkar  and  Per la  Subbaiah 
University  of  Rochester  and  Oakland  University 


. 


* Research  sponsored  by  the  Air  Force  Office  of  Scientific  Research,  Air 
Force  Systems  Command,  USAF  under  Grant  No.  AFOSR-77-3360.  The  United 
States  Government  is  authorized  to  reproduce  and  distribute  reprints 
for  Governmental  purposes  notwithstanding  any  copyright  notation  hereon. 
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TABLE  6 


EXACT  POWER  FUNCTioN  OF  HOTEl LING’S  T"  TEST  PROCEDURE* 
a = .01  and  .05,  p = 2 


1.5 

1.00 

1.00 

1.00 

1.00 

1.2 

.995 

.996 

.996 

.999 

1.0 

.965 

.966 

.970 

.989 

0.8 

.845 

.851 

.867 

.943 

0.5 

. 43b 

.451 

.495 

.742 

0.2 

. 105 

.119 

. 165 

— 

0. 1 

.063 

.076 

.466 

0.0 

.050 

.050 

.237 



.018 

.032 

.206 

.010 

.014 

.027 

.195 

0.0 

0.1 

0.2 

0.5 

1.00 

1.00 

1.00 

1.00_^ 

' 

1.00 

1.00 

1.00^ 

......^ — ' — 

.998 

1.00 

1.00 

1 . 5 

.990 

^ 

1.00 

1.00 

1.2 

' 

.995 

.999 

1.00 

1.0 

.935 

.931 

.996 

1.00 

0.8 

.790 

.929 

.984 

.999 

0.5 

.640 

.863 

.967 

.998 

0.2 

.614 

.851 

.963 

.998 

0.1 

.605 

.846 

.962 

.998 

0.0 

0.8 

1.0 

1.2 

l.S 

I. 


I 


u 


1 


TABLE  7 


POWER  FUNCTION  OF  FISHER'S  COMBINATION  OF  THE  STEP-DOWN 
TESTS  ESTIMATED  FROM  THE  MONTE  CARLO  EXPERIMENT 

a = .01  and  p = 2 


1 . 3 

.997 

.998 

.998 

.999 

1.00 

1.00 

1.00 

1.00 

1.2 

.958 

.950 

.963 

.990 

.997 

1.00 

1.00 

1.00 

1.0 

.815 

.324 

.837 

.928 

.983 

.996 

1.00 

1.00 

0.3 

. 577 

.582 

.623 

.794 

.949 

.981 

.998 

1.00 

0.5 

.177 

. 194 

.221 

.496 

.300 

.928 

.983 

.999 

0.2 

.023 

.030 

.049 

.241 

.655 

.350 

.963 

.998 

0. 1 

.014 

.023 

.034 

. 198 

.593 

.839 

.953 

.996 

0.0 

.011 

.013 

.028 

.185 

.594 

.844 

.962 

.997 

0.0 

0.1 

0.2 

0.5 

0.8 

1.0 

1.2 

1 . 5 

u 

1 

TABLE  8 


POWER  FUNCTION  OF  TIPPETT'S  COMBINATION  OF  THE  STEP-DOWN 


TEST 

ESTIMATED 

FROM 

THE  MONTE  CARLO 

EXPERIMENT 

a 

= .01 

and 

p = 2 

1.5 

.998 

.997 

.998 

.995 

.997 

.998 

1.00 

1.00 

1.2 

.965 

.954 

.956 

.954 

.962 

.990 

.998 

.999 

1.0 

.863 

.860 

.839 

.835 

.916 

.965 

.995 

.999 

0.8 

.614 

.614 

.615 

.627 

.813 

.933 

.989 

1.00 

0.5 

.203 

.209 

.201 

.349 

.691 

.893 

.977 

1.00 

0.2 

.025 

.027 

.041 

.226 

.676 

.876 

.978 

.999 

0.1 

.017 

.019 

.030 

.204 

.644 

.883 

.973 

.998 

0.0 

.011 

.010 

.030 

.208 

.656 

.885 

.981 

.999 

0.0 

0.1 

0.2 

0.5 

u 0-8 
1 

1.0 

1.2 

1.5 

Below-diagonal  elements  correspond  to  a * .01  and  above-diagonal 
elements  correspond  to  a » .05 
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TABLE  9 

POWER  FUNCTIONS  OF  T‘‘-TEST  AND 
MODIFIED  STEP-DOWN  TESTS 
a = . 01  , p = 3 


Configuration 


0.0 

0.1 

0.2 

0.4 

0.6 

0.8 

1.0 

Fisher 

.013 

.017 

.066 

.366 

.818 

.987 

.999 

Tippett 

.010 

.016 

.049 

. 198 

.518 

.826 

.972 

T" 

.012 

.018 

.062 

.340 

.791 

.981 

.999 

T" -exact 

.010 

.019 

.058 

.344 

.795 

.980 

.999 

Fisher 

.010 

.012 

.025 

.085 

.243 

.498 

.753 

Tippett 

.009 

.009 

.025 

.092 

.298 

.600 

.846 

T^ 

.011 

.021 

.024 

.080 

.237 

.498 

.758 

T^-exact 

.010 

.013 

.022 

.080 

.235 

.491 

.753 

TABLE  10 

POWER  FUNCTIONS  OF  T^-TEST  AND 
MODIFIED  STEP-DOWN  TESTS 
ct  = .01  , p = 4 


Configuration 

Tests 

u 

0.0 

0.1 

0.2 

0 .4 

0.6 

0.8 

1.0 

u 

Fisher 

.012 

.027 

.066 

. 458 

.386 

.996 

1.00 

Tippett 

.012 

.018 

.043 

.195 

.491 

.808 

.966 

Iv 

T” 

.010 

.027 

.065 

. 422 

.854 

.993 

1.00 

U 

"7 

T"-exact 

.010 

.019 

.063 

.400 

.860 

. 992 

1.00 

Fisher 

.015 

.017 

.030 

.072 

.196 

.400 

.634 

1 A 

Tippett 

.012 

.012 

.021 

.081 

.262 

.536 

.816 

\ V 

T“ 

.015 

.014 

.026 

.060 

.184 

.389 

. 636 

0 

T^-exact 

.010 

.012 

.019 

.063 

.182 

.400 

.659 
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